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ABSTRACT 

We analyze six X-ray bright elliptical galaxies, observed with Chandra and XMM- 
Newton, and approximate their gravitational potentials by isothermal spheres ip = 
T;^logr over a range of radii from '^0.5 to ~ 25 kpc. We then compare the circular 
speed Vc derived from X-ray data with the estimators available from optical data. In 
particular we discuss two simple and robust procedures for evaluating the circular 
speed of the galaxy using the observed optical surface brightness and the line-of- 
sight velocity dispersion profiles. The best fitting relation between the circular speeds 
derived from optical observations of stars and X-ray observations of hot gas is Vc,opt — 
V Vc,x, where rj = 1.10 — 1.15 (depending on the method), suggesting, albeit with 
large statistical and systematic uncertainties, that non-thermal pressure on average 
contributes ~20-30% of the gas thermal pressure. 

Key words: Galaxies: Kinematics and Dynamics, X-Rays: Galaxies: Clusters 



1 INTRODUCTION 

In spiral galaxies, disk rotation curves offer an accurate 
and robust way of measuring total gravitational potentials 
to distances as large as 10-30 kpc. To a first approxi- 
mation the rotation curves are flat over a broad range 
of radii, suggesting an isothermal (logarithmic) potential 
characterized by ip{r) = v^logr. In early-type galaxies, 
measuring the gravitational potential is much more difficult 
since there are no tracers such as cold gas or disk stars on 
orbits of known shape. Several methods have been used 
to measure the potentials of elli ptical galaxies includin; 



ies (iKooDmans et al. 


l2006l; iMandelbaum et all 20061: 


Gavazzi et al.l 2007; 


Mandelbaum, van de Ven, & KeetonI 


2008]) and ffor the most massive galaxies) modeling the 



hydr ost atic atmospheres of X-ray emitt in g gas ([Mathews 



19781 : iForman, Jones fc Tucker! 1 19851: iFukazawa et al 



20od : iHumphrev et al.H2006l : IChurazov et aLlbOOSl ). Recent 



20091 ). 



2001 



detailed modeling of s t ellar orbits jKronawitter et al.. 

Cappellari et al.l l2006l : iThomas et al.l l2007l : iGiiltekin et al I 



studies using stellar kinematics and lensing suggest that the 
pote ntials of early-type ga l axies are approximately isother- 
mal jGerhard et al.H200ll : iTreu et al.f l2006l : iGavazzi et al] 
|2003), similar to disk galaxies. Consistency of the mass 
profile with M(r) oc r (implying an isothermal potential) 
was also suggested for several ellipt ical galaxies based on 



tracers such as globular 
and 



clusters, planetary neb- 



ulae, and satellite galaxies ( Romanowsky fc Kochanekl 
I2OOII : ICoccato et al] l2009l : iKlvpin fc Prada ,2009, ) . strong 



19851: Buotc & Canizared 


19941: Kim & Fabbiand 


1995: 


Nulsen & Bohrinac3 1995 


: Buote & Canizaresl 


1998: 


Fukazawa et al.l bood). Here we present independent evi- 
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dence that the potentials of bright elliptical galaxies are 
close to isothermal, from X-ray observations with Chandra 
and XMM-Newton. If the potential of a typical elliptical 
galaxy is indeed not far from being isothermal (for a 
range of radii) then it can be characterized in that radius 
range with a single number - the circular speed Vc. This 
makes the comparison of potentials derived from X-ray 
and optical data especially simple since it does not require 
a point-by-point comparison and for a pure isothermal 
(logarithmic) potential the results are not sensitive to the 
range of radii used for evaluation of Vc- While deviations 
from isothermality are certainly present at some level, we 
believe this approach is useful. In this paper we develop 
simplified methods for characterizing the X-ray and optical 
data and comparing the results to place constraints on the 
non-thermal pressure in the hot gas in elliptical galaxies. 

By construction our method is quick and approximate 
and is not intended to replace a careful and comprehensive 
analysis of individual objects. It might be useful, for exam- 
ple, in a larger sample, when detailed modeling is not prac- 
tical due to noisy or missing data. We illustrate the method 
on a small and rather arbitrarily selected sample of X-ray 
bright elliptical galaxies. 

The structure of this paper is as follows: in iJ5] we de- 
scribe our sample of galaxies and how we derive the gravita- 
tional potential from X-ray observations, and in !j3]we dis- 
cuss methods for determining the potential from optical ob- 
servations of stellar velocity dispersions. The implications of 
our results for the mass distribution and non-thermal pres- 
sure contribution in these galaxies are discussed in and 
^contains conclusions. 



2 THE SAMPLE AND X-RAY ANALYSIS 

For our analysis we selected six nearby (distance less than 
~30 Mpc, see Table[T]), X-ray bright galaxies, which were all 
well observed with Chandra and XMM-Newton. All galaxies 
in the sample are very bright and dominate (in terms of 
mass or potential) their environment up to at least several 
effective radii. This (at least partly) justifies the standard 
practice of deriving mass/potential profiles from X-ray data 
using the assumption that the X-ray emitting gas forms a 
hydrostatic atmosphere. 

For the analysis we used publicly available Chandra and 
XMM-Newton data. Combining the data from these two in- 
struments provides a cross-check of the results and also gives 
better constraints on the innermost and outermost regions, 
thanks to the superb angular resolution of Chandra and the 
large field of view of XMM-Newton, respectively. 

For Chandra t he data were prepared following the pro- 
cedure described in lVikhlinin et ahl (|2005l V This includes fil- 
tering of high background periods and application of the 
latest calibration corrections to the detected X-ray photons, 
and determination of the background intensity in each ob- 
servation. 

For XMM-Newton the data were prepared by remov- 
ing background flares using the light curve of the detected 
events above 10 keV and re-normalizing the "blank fields" 
background to match the observed count rate in the 11-12 



keV band. In the subsequent analysis we use the data from 
the EPIC/MOS detector only. 

The analysis of the X-ray data is based on 
a non-paramet r ic deprojec tion procedure, described in 
I Churazov et all (|2003l . |2008| , COS hereafter). In brief, the 
observed X-ray spectra in concentric annuli are modeled 
as a linear combination of spectra in spherical shells; the 
two sequences of spectra are related by a matrix describing 
the projection of the shells into annuli. To account for the 
projected contribution of the emission from shells at large 
distances from the center (i.e., at distances larger than the 
radial size Vmax of the region well covered by actual ob- 
servations) one has to make an explicit assumption about 
the behavior of the gas density/temperature profile at large 
radii. We assume that at all energies the gas volume emis- 
sivity at radii beyond Tmax declines as a power law with 
radius. The slope of this power law is estimated based on 
the observed surface-brightness profile within the range of 
radii covered by observational data. Since we assume that 
the same power law shape is applicable to all energy bands, 
effectively this assumption implies constant spectral shape 
and therefore the isothermality of the gas outside rmax (here 
"isothermal" means that for r ^ rmax the gas temperature 
is independent of radius, not that the gravitational poten- 
tial is logarithmic). The contribution of these layers is added 
to the projection matrix with the normalization as an addi- 
tional free parameter. While the limitations of this approach 
are obvious, the contribution of these outer shells is usually 
important only in the few outermost radial bins inside Tmax, 
especially when the surface-brightness profile is steep. The 
final projection matrix is inverted and the shells' spectra are 
explicitly calculated by applying this inverted matrix to the 
data in narrow energy channels. 

The resul ting spectra are approximated in XSPEC 
I ArnaudI [l996l ) with the AP EC one-temperatur e optically 
thin plasma emission model jSmith et al.ll200ll ). The red- 
shift z (from the NASA/IPAC Extragalactic Database - 
NED) and the line-of-sight column density of neutral hy- 
drogen Nh (based on Kalberla et al. 2005) have been fixed 
at the values given in Table [T] For each shell we determine 
the emission measure (and therefore gas density) and the 
gas temperature. These quantities are needed to evaluate 
the gravitational potential through the hydrostatic equilib- 
rium equation. For cool (sub-keV) temperatures and approx- 
imately solar abundance of heavy elements, line emission 
provides a substantial fraction of the 0.5-2 keV fiux. With 
the Chandra and XMM-Newton spectral resolution the con- 
tributions of continuum and lines are difficult to disentangle. 
As a result the emission measure and abundance are anti- 
correlated, which can lead to large scatter in the best-fit 
emission measures. As an interim (not entirely satisfactory) 
solution, we fix the abundance at 0.5 sola r for all shells, using 
the de fault XSPEC abundance table of I Anders fc Grevessd 
(|l989l ). We return to this issue in 

With known gas density n and temperature T in each 
shell we can use the hydrostatic equilibrium equation to eval- 
uate the gravitational potential ip: 

p dr dr 
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Table 1. Sample of elliptical galaxies. The columns are: (1) - common name of the galaxy; (2) - re dshift from the NASA/IPAC 
Extragalactic Database; (3) - adopted distance; (4) - hydrogen column density from lKalberla et al.l | |2005| '); (5) effective radius (eg. I14II . 
(6) Sersic index, (7) - central velocity dispersion standardized to ~ 0.6 kpc aperture from Hyperleda, except for NGC 4472, taken from 
B94, (8) - line-of-sight velocity dispersion at the "sweet spot" Rs; (9) - circular speed estimated from the central velocity dispersion 
as fc.c = \/2o"c, (10) - circular speed Wc,s accordi ng to i)3.4.1l (11) - circular speed v^- i accordin g to il3^.2l References for the effectiv e 



radius, Sersic index and ste llar kine matics: [D941 -iD'O nofrio. Capacci oli. &: CaonI |l9 94[): [C93] -icaon. 

[SOS] - ISpolaor et~all bOOSi): [M05] -iMahdavi. Trentham. fc TuUvl bOOSl ) ; [KO O] -iKronawitter et al.l 1I2OOOI); fSOOl ^ 
fB94] - iBender. Saglia. fc Gerhard! lll994): fG 09] - iGebhardt fc Thomad l|2009l 'l: iMurhphv fc Gebhardtl l|2009t l. [DDI] 



Capaccioli. fc D'Onofrio 
ISagU^'et'af 



1I2OOII ). Distances are from lTonrv et al.l fooib . Note that distances are not explicitly used in the subsequent analysis 



1 



1993 ); 



200c ); 



De Bruvne et al.l 



Name 


z 


D, Mpc 


Nh 




Re, arcsec 


Sersic index, n 


CTc 


<j{Rs) 


Vc,c 


Vc,s 




(1) 


(2) 


(3) 


(4) 




(5) 


(6) 




(7) 


(8) 


(9) 


(10) 


(11) 


NGC1399 


0.00475 


20.0 


1.5 


1020 


117 


12.24 


[D94] 


341 


242 [SOO] 


482 


412 


394 


NGC1407 


0.00593 


28.8 


5.4 


1020 


70 


8.35 


[S08] 


272 


256 [S08] 


385 


435 


408 


NGC4472 (M49) 


0.00333 


16.3 


1.5 


1020 


257 


6.27 


[C93] 


320 


289 [B94] 


452 


492 


445 


NGC4486 (M87) 


0.00436 


16.1 


1.9 


1020 


145 


6.51 


[C93] 


336 


312 [G09] 


475 


530 


536 


NGC4649 (M60) 


0.00373 


16.8 


2.0 


1020 


118 


5.84 


[C93] 


336 


244 [DOl] 


475 


414 


436 


NGC5846 


0.00572 


24.9 


4.3 


1020 


79 


3.95 


[M05] 


241 


215 [KOO] 


341 


366 


338 



Table 2. Best fitting Vc,x derived from Chandra and XMM-Newton data on the potential profiles approximated with a logarithmic law 
ip(r) = j;2 logr + b. The quoted uncertainties are pure statistical errors, determined from a Monte Carlo procedure. The last column is 
the average of the results from Chandra and XMM-Newton in the preceding two columns. 



Name 


ri{') 


r2{') 


f'c.Chandrakm S ^ 


VcXMM km s ^ 


Vc^x km s ^ 


NGC1399 


0.1 


5.0 


403 ± 1.8 


395 ± 2.1 


399 


NGC1407 


0.1 


2.0 


368 ± 9.8 


356 ± 9.8 


362 


NGC4472 (M49) 


0.1 


5.0 


372 ± 4.0 


367 ± 2.0 


370 


NGC4486 (M87) 


0.1 


5.0 


448 ± 1.3 


437 ± 1.6 


443 


NGC4649 (M60) 


0.1 


5.0 


417 ± 2.4 


422 ± 1.8 


420 


NGC5846 


0.1 


5.0 


335 ± 3.1 


331 ± 5.2 


333 



where p — ^nipn is the gas density, P = nkT is the pres- 
sure, /i is the mean atomic weight of the gas (p = 0.61 
assumed throughout the paper), irip is the proton mass and 
k is the Boltzmann constant. Integrating the above equa- 
tion one gets an expression for the gravitational potential 
through the observables n and tQ 



T 



dlogn 
dr 



dr + T 



+ C, 



(2) 



where C is an arbitrary constant. We choose the constant C 
such that If [Re) ~ 0, where Re is the optical effective radius 
(see eq.[T3]and Table [T]). The resulting potentials are shown 
in Figure [1] with red (Chandra) and XMM-Newton (blue) 
points. The error bars were evaluated by a Monte Carlo 
procedure, starting from the measured values of n and T in 
each shell, adding "noise" to the data points and re-deriving 
the potential via equation ^ (see COS for a discussion of 
the limitations of this procedure). 

As discussed in COS the expression ^ for the potential 
can be evaluated directly from deprojected data - density 
and temperature profiles. Since the logarithmic derivative 
of the density is under the integral the enhancement 

of errors due to differentiation is not a major issue and 
evaluation of eq. (HJ does not require any parametric de- 



scription of the density and/or temperature profiles. This 
removes the ambiguity in the choice of a functional form 
for these profiles. Since the potential is calculated directly 
via eq. Q no apriori parametrization of the potential (or 
mass) is required. If information on the gravitational poten- 
tial is available from other data (e.g. from optical observa- 
tions) the comparison of the potential profiles can be done 
directly with the results of eq. ([2|. In COS the relation be- 
tween potential profiles derived from X-ray and optical data 
was written as ipx = afopt + b and the value of a was evalu- 
ated. In this case any transformation of the potential profile 
to the circular speed (or equivalently to the mass profile) 
which involves differentiation of the potential profile is not 
needed and should be avoided. 

Thus, the application of eq. ((2)1 to the X-ray data is 
essentially free from any parametrization and does not rely 
on direct differentiation of the data. This is the main dif- 
ference with most of the other techniques which either use 
parametrization of the density /temperature profiles or make 
an assumption on the form of the mass profile (see e.g. 



Mathewslll978l:lF abian ct al."l98ll: lForman. Jones, fc Tucker] 



^ Throughout this paper, "log" denotes natural logarithm. 



19851: iNialsen fc B ohrinacr 1993: iHumphrev et al.1 l2006l : 

Fukazawa et al Eood for various techniques of the mass pro- 
fil e reconstruction). The m ost close approach to ours is that 
of [Humphrey et al] (|2006l ). where the temperature and mass 
profiles are parametrized, the hydrostatic equilibrium equa- 
tion is solved to find the density and the best-fitting param- 
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Figure 1. Gravitational potentials derived from Chandra (red) and XMM-Newton/MOS (blue) observations of each galaxy. Potentials 
are normalized to zero at i?e . Vertical lines mark the range of radii used to approximate the data by a log r + b law, and the thin solid 
lines are the best-fit approximations of the Chandra and XMM-Newton data to this law. The vertical dotted lines mark the effective 
radius. 



eters are iteratively found. As mentioned above in eq. ((Sjl all 
observables are on the r.h.s and one can reconstruct the po- 
tential directly and defer the parametrization (if needed) to 
a final step of manipulations with reconstructed potential. 

We are now considering the case when limited infor- 
mation on the potential from optical data is available. For 



instance, consider a situation when only central velocity dis- 
persion of the galaxy is known (see discussion in In this 
case to compare X-ray and optical potentials a parametriza- 
tion of the potential derived from eq. ([2]) is needed. An es- 
pecially simple and straightforward parametrization is pos- 
sible if the potential is isothermal (p = v1 log r -\- b, since in 
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this case a single number ~ circular speed Vc can character- 
ize the potential. In fact, recent studies using stellar kine- 
matics and lensing do suggest that the pote ntials of early- 
type galaxies are approximately isotherma l (|Gerhard et alj 
l200ll : lTreu et all 120061 : iGavazzi et al.ll2007l '). similar to disk 
galaxies. Consistency of the mass profile with M(r) cx r 
(implying an isothermal potential) was also suggested for 
several elli ptical galaxies based on the analysis of X-Ray 
data ( e.g.. iTrinchieri fc Fabbiandll985l: iBuote fc Canizared 



1994 iKim fc Fabbiand I1995I: iNulsen fc BohringeJ Il995l : 



Buote &: Canizared 19981 : Fukazawa et al.1 20061 ). Below we 
confirm an approximate isothermality of the potentials us- 
ing the results of application of the non-parametric method 
(eq. [2]) and determine the best-fitting value of Vc- 

The axes in Figure [T] are log-linear and therefore an 
isothermal (logarithmic) potential should look like a straight 
line log r + b. To first order this is true, although there 
are statistically significant deviations, particularly at the in- 
nermost and outermost radii, which we discuss below. The 
agreement of Chandra and XMM-Newton data is good, ex- 
cept for the inner region where the better spatial resolution 
of Chandra is important. We choose to ignore the data in- 
side the ri = 0.1' circle where this effect is apparent. We 
also introduce a cutoff at large radii (typically r2 ~ 5'), 
where the results are sensitive to the assumed extrapola- 
tion of the emissivity profile. The actual values of ri and 
r2 used in the analysis are given in Table [2] and are shown 
in Figure [T] as a pair of vertical lines. Between ri and r2 
the potential was approximated with a logarithmic function 
log r + b with Vc and b being free parameters. Best fitting 
values were found by minimizing the root-mean-square de- 
viation between the model and observed potential profiles. 
The resulting values of Vc,x for Chandra and XMM-Newton 
data are given in Table [2] as Dc, Chandra and Wc.xmm- The 
statistical uncertainties determined from the Monte Carlo 
procedure are also given in Table [S] In general there is 
good agreement between the values obtained by the two in- 
struments, confirming that the uncertainties introduced by 
statistical errors and cross-calibration uncertainties between 
the two instruments (including background subtraction pro- 
cedures) are small. One can expect however that the real un- 
certainties are dominated by systematic errors arising from 
our model assumptions (e.g., the assumption of spherical 
symmetry used in the deprojection analysis), which affect 
both datasets in similar ways. For subsequent analysis we 
use the average of the results from the two instruments, 

Vc,X = {Vc jChandra 

+ Uc,xmm)/2 (last column in Table[2l). 

The deviations of the potential from the isothermal 
shape (f = alogr + b can be studied by assuming a func- 
tional form y = ar^ + b and looking for the best-fitting value 
of a. For the five galaxies with ri — 0.1' and r2 = 5' the 
best-fitting value of a varies from 0.03 for NGC1399 to 0.16 
for NGC4472 (the mean value of a is 0.11). For NGC1407 
(ri — 0.1' and r2 = 2') the value of a is ~ 0.5. This sug- 
gests that on average the profiles are slightly concave, i.e. 
circular speed slightly increases at large radii, which can be 
seen in Fig. [T] This is consistent with the fact that these 
massive ellipticals are sitting at the centers of more massive 
group/cluster size halos which dominate at large radii. Given 



the mean value of a = 0.11 one can estimate that changing 
both ri and 7-2 by a factor of two would on average change 
the estimate of Vc by a factor of only 2°'^'^ ~ 4%. 

We explicitly tested the sensitivity of Vc to the values of 
ri and r2 for all objects in the sample by increasing the lower 
boundary by a factor of two (i.e. ri = 0.2') and recalculating 
Vc- The results of this test are given in Table [3] (column 
Arix2), where 



Vc,ii x2 — Vc,X 
Vc,X 



(3) 



and Vc,x is the circular speed from Table [S] and Wc,rix2 is 
calculated as a mean value of circular speeds measured using 
Chandra and XMM-Newton for increased ri (similarly to 
Vc,x in Table [2]). In a separated test we decreased the value 
of r2 by a factor of two (i.e. r2 = 2.5' for all objects except 
for NGC1407, where r2 — 1') and again recalculated Vc (see 
TableO column A^^^2)- From Table[3]it follows that factor 
of 2 changes in either ri or r2 causes few per cent changes 
in Vc. 



2.1 Flat abundance profile 

We now illustrate the impact of our assumption of a flat 
abundance profile, using NGC1399 as an example. 

As is obvious from equation ((5}, the absolute normal- 
ization of the gas density n does not affect the calculations 
of the potential. From this point of view it is not the partic- 
ular value of the heavy-element abundance in the spectral 
models, but rather the radial variation of the abundance 
that affects the derived potential profile. In elliptical galax- 
ies one can expect an increase of metal abundance towards 
the center of the galaxy. This is usually true for a range of 
radii except for the very center, where a "dip" in the metal 
abundance is often seen when fitting the data with a single 
temp erature plasma emission model (e.g.. [Matsushita et al.l 
120021) . 

As mentioned above, measuring the metal abundance 
from X-ray spectra in cool systems is difficult because of 
the ambiguity of separating line from continuum emission 
with the limited energy resolution of X-ray CCDs. For a 
multi-temperature plasma this separation is even more com- 
plicated because fitting the emission with one-temperature 
models leads to a biased estimate of the abundance (e.g., 
lBuotd[2000l ). In deprojected spectra, such as we use here, 
the bias arising from a multi-temperature plasma is reduced 
because the superposition of cooler and hotter emission com- 
ing from different radii is removed (provided the object is 
spherically symmetric), although if the plasma is intrinsi- 
cally multi-temperature then the problem remains. On the 
other hand the signal-to-noise of the deprojected spectra is 
much lower than for the projected spectra. It is therefore 
desirable to keep the number of free parameters in the fit as 
small as possible. 

As an example we show in Figure [2] (left) the gas param- 
eters (electron density, temperature and abundance of heavy 
metals) as a function of radius for NGC1399. In two models 
(blue circles and green triangles) the abundance was fixed at 
0.5 and 1 times solar, respectively. In the third model (red 
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squares) the abundance was a free parameter. One can see 
that in the third model the errors on the abundance are sub- 
stantial and, as expected, the abundance and the gas density 
are anti-correlated. The potential profiles corresponding to 
these spectral models are shown in the right panel of Figure 
[2] Clearly there is a significant change in the potential curves 
(curves are normalized to zero potential at 1'.5). Formally 
calculated values of Vc for the three spectral models are: 



Vc = 421 km s" 
Vc = 403 km s"^ 
= 396 km s"^ 



V, 



free abundance 
abundance = 0.5 x solar 
abundance = 1 x solar. 



(4) 



Thus, there is a substantial, but not dramatic effect of the as- 
sumed abundance profile on the derived logarithmic slope of 
the potential profile. Given that the limited statistics in our 
deprojected spectra do not allow for a robust abundance de- 
termination for all galaxies in the sample, we decided to keep 
the assumption of a flat abundance profile so that we could 
analyze all objects in a uniform way. It is likely that this 
approximation introduces errors in the best-fitting value of 
Vc of roughly 5-15 km s~^. If all galaxies in our sample have 
a similar abundance profile to the one derived in NGC1399, 
one can expect the values of Vc obtained under the assump- 
tion of a fiat profile with 0.5 x solar metallicity to be biased 
low by 3-4%. 

Relative changes in the circular speed when abundance 
of heavy elements is a free parameter (constrained to be 
in the range from to 2) for all objects in the sample are 
calculated in Table [3] (column Aabund): 



^c.frcc abunda 



- Vc,X 



Vc,X 



(5) 



where Vc,free abundance IS Calculated as a mean value of cir- 
cular speeds measured using Chandra and XMM-Newton, 
similarly to Vc,x (see ij2|. 

Yet another parameter related to the chemical com- 
position of the hot gas that is important for the eval- 
uation of the gravitational potential is the mean atomic 
weight jj, (see eg. [H. For the so l ar ph otospheric abun- 
dance table of lAnders fc Grevessj (| 19891 ') fi = 0.614 (as- 
suming fully ionized plasma). If we lower the abundance 
of heavy elements to 0.5 solar, fi changes to 0.610, i.e. by 
0.7%. Changing the helium abundance would, of course, 
have a much stronger impact on /i. For instance, doubling 
the abundance of helium relative to the solar photospheric 
mix would increase jj, to 0.695, i.e., by 13%. One mech- 
anism that could lead to helium enrichment beyond so- 
lar abundan ce is gravitational sedimentation of heavy el- 
ements (e.g.. [ Cilfanov fc Sunvaevlll984l: IChuzhov fc NusseJ 



l2003l : IChuzhov fc Loebll2004l : lEttori fc FabianlbOOd V If sed- 
imentation is indeed important for the interstellar gas in the 
cores of elliptical galaxies, then it will also affect the emis- 
sivity. In the present work we assume that jj. = const — 0.61 
at all radii, but one should bear in mind that the results 
might change if the helium abundance changes with radius 
or simply differs from the canonical value of 7.92 x 10~~ 
(number density relative to hydrogen) adopted here. 



2.2 Deviations from spherical symmetry 

To make a crude estimate of an impact on Vc by possible 
deviations from spherical symmetry one can make deprojec- 
tion analysis in individual wedges and compare the results. 
We divided the data on each object into two independent 
halves - Northern and Southern and repeated the analysis 
for each half separately. The difference in the circular speeds 
between two halves is given in Table [3] (column Ans ) : 



Vc.N — Vc,S 
Vc,X 



(6) 



where Vc,n and Vc,s are the circular speeds for the Northern 
and Southern halves respectively. For most of the objects 
the difference is at the level of 2-4% (compared to ~1% 
pure statistical uncertainty). 

The largest difference ~8% is for NGC1407. Inspection 
of the data have shown that this difference primarily comes 
from XMM-Newton data. We note that NGC1407 was ob- 
served by XMM-Newton with an offset angle of ~ 8' and 
the PSF distortions might contribute to this difference. 

In the analysis of optical data (|3]) we discuss optical 
constraints on the gravitational potential, assuming that the 
galaxies are spherically symmetric. All objects in the sample 
are round, being E0-E2 galaxies, although in some objects 
(e.g. M87) the ellipticity increases with radius. We note here 
that in the radial range where stars are dominating the mass, 
the potential is more spherical than the distribution of the 
stars themselves. The same is likely true even in outer re- 
gions, given a steep decrease [R^^ or steeper) of the optical 
surface brightness. 

The analysis of independent wedges is not a useful indi- 
cator of a possible error in the effective circular speed when 
the object posses an axial symmetry and is viewed pole-on. 
Possible effects of a non-spherical potential on optical es- 
timates for a pole-on galaxy are discussed in C08 (Section 
7.2 there). There we show that for an axis ratios (along 
and perpendicular to the line-of-sight) of 0.5 (oblate) and 2 
(prolate) the error in mass is factor of 0.79 and 1.04 respec- 
tivelj|3, and argue that the probability that many objects 
in our sample are strongly oblate or oblate systems viewed 
pole-on is rather low. Therefore on average the error is going 
to be smaller than these extreme values. 

X-ray analysis of the ellipsoidal objects is discussed 
in iPiffaretti. Jetzer. fc Schindled |2003| and C08. If poten- 
tial is logarithmic and the isothermal gas is in hydro- 
static equilibrium, then the potential is recovered cor- 
rectly from the spherically symmetric analysis (see Ap- 
pendix B in C08) . For a more complicated potential or non- 
isothermal gas the deviations are p resent (see e.g. Fig.l in 
IPiffaretti. Jetzer. fc Schindleij uOOSl . for an example of ex- 
pected bias for prolate and oblate cases in A2390 cluster). 
Although the effects of non-spherical potential are differ- 
ent for X-ray and optical analysis, the sign of the effect for 
oblate and prolate systems is the same. Therefore errors in 
mass/potential partly compensate each other when effective 
circular speeds are compared. 



See COS for the assumptions used in this calculation 
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Figure 2. Left: Radial profiles of deprojccted gas parameters (abundance, temperature, electron density) in NGC1399. The parameters 
were obtained using a single-temperature APEC model (with fixed low-energy absorption and redshift) fit to the deprojected Chandra 
spectra from a set of spherical shells. For the blue and green models (circles and triangles) the abundance was fixed at 0.5 and 1 
times solar, respectively, while for the red model (squares) the abundance was a free function of the radius. Right: Potential profiles 
corresponding to the spectral models shown in the left panel. Straight lines show best-fitting approximations ip = v'^ log r + b. 



We concluded that the assumption of a spherical po- 
tential can contribute to the discrepancy between X-ray and 
optical data, but it is unlikely that on average the magnitude 
of the effect exceeds several per cent. 



2.3 Contribution of unresolved LMXBs and weak 
sources 

The contribution of low-mass X-ray binaries (LMXB) can 
dominat e the X-ray emission in gas -poor elliptical galax- 
ies (e.g., iTrinchieri fc Fabbiand[l985^ . Even when the spa- 
tial resolution and sensitivity allows the resolution of some 
LMXBs the remaining emission can still be dominated by 
even weaker unr esolved sources, such as accreting white 
dwarfs (see, e.g.. iRevnivtsev et al.|[2008l . for the analysis of 
the galaxy NGC3379). 

Since all objects in our sample are gas-rich galaxies the 
contribution of unresolved sources is not a major issue, espe- 
cially for such bright objects as M87 and NGC1399. For less 
luminous objects, because of their harder spectrum com- 
pared to the hot gas emission, LMXBs can affect the ob- 
served fluxes at the high end of the Chandra and XMM- 
Newton energy bands. To test the magnitude of this ef- 
fect we added a power law component with a photon index 
r = 1.6 with a free normalizatio n to represent the contri- 
bution of unresolve d source s (e.g., Irwin . Athev. fc BregmarJ 
I2OO3I : iDavid et al.1 2006; H umphrey fc BuotQ,2006i) and re- 



calculated the circular speed for NGC4472. The resulting 
values are only ~1% lower than in Table [2] 

Relative changes in the circular speed when a power 
law component with a photon index T = 1.6 is added to the 
spectral model for all objects in the sample are calculated 
in Table 13] (column Almxb): 



t'c.LMXB — Vc,X 
Vc,X 



(7) 



On average adding a power law component shifts the circular 
speed few percent lower. The largest effect is for M87 and in 
this very gas rich galaxy the change in Vc likely reflects the 
complexity of the spectrum rather then the contribution of 
real LMXBs. 



2.4 Summary on uncertainties in X-ray analysis 

The summary of the uncertainties in determination of the 
circular speed from X-ray data is given in Table [3] The 
columns in the table show the relative deviation of the circu- 
lar speed from the reference value Vc,x given in Table[5]when 
changes are made to the analysis procedure. The quoted un- 
certainties in columns labeled Err. are pure statistical errors. 
They have clear meaning only when the difference between 
the Northern and Southern parts of the galaxies are con- 
sidered, since the data are independent. In all other cases 
the quoted uncertainty corresponds to the largest statistical 
error in one of the two values of Vc used to calculate A. 

The last two rows in the Table[3]give the mean change in 
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the circular speed and the RMS value (relative to zero). The 
values given in the Table allows one to get an idea of the 
uncertainties introduced by the assumptions incorporated 
into calculations of Vc. For instance, letting abundance be a 
free parameter on average shifts the value of Vc up by ~ 2%, 
while adding a power law component (to control possible 
contribution of LMXBs) on the contrary shifts Vc down by 
~ 2%. Changing any of the bounds of the fitting range ri 
and r2 by a factor of 2 changes Vc by few % up or down. 

While rigorous evaluation of combined uncertainties is 
difficult, we can very crudely estimate it by adding RMS 
values of Aabund, Almxb, Ans, Ari and Ar2 quadratically. 
The resulting value 7.1% is crude a characteristic of the un- 
certainties in the final value of Vc,x introduced by the modi- 
fications of our analysis procedure. Note that pure statistical 
errors also contribute to the above estimate. 



3 ROBUST ESTIMATORS OF THE 
ISOTHERMAL POTENTIAL FROM 
OPTICAL DATA 

The results of the X-ray analysis suggest that the shape 
of the potential profile of the galaxies in our sample is not 
far from isothermal over the range of radii ~ 0.5-25 kpc, 
corresponding to ~ 0.05-3 Re- Let us assume that this is 
true at all radii and the total gravitational potential of a 
galaxy is logarithmic. 



(y3(r) = log r + const. 



(8) 



It is natural to ask, given this assumption, how the velocity 
dispersion of the stars in the galaxy is expected to be re- 
lated to Vc- A robust method to determine Vc from optical 
observations would enable us to check our results by com- 
paring the derived values in Table [J] to determine Vc and 
thus the gravitational potential for many galaxies without 
X-ray observations, and to look for evidence of non-thermal 
pressure support in the X-ray gas. 

For simplicity we shall assume that the galaxy is spher- 
ical with known surface brightness I{R) and line-of-sight 
velocity dispersion (j{R), and ask what is the best way to 
determine Vc from known I{R), cj{R) given that the poten- 
tial is isothermal. 



3.1 Circular speed from velocity dispersion for 
distant (unresolved) galaxies 

Since the acceleration of a star in the logarithmic potential 
i' = — — v^r/r'^, the virial theorem for such a system 
reads 



= (r^ -f- r • r) = (v^) — 



(9) 



where v = r. Since the galaxy is spherical, the mean-square 
velocity along any line of sight is equal to |(v'^) and there- 
fore 



^^irdRRl{RW{R) 



(10) 



J^dRRI{R) ' 

i.e., the circular speed is equal to ^3 times the integrated 
line-of-sight dispersion for the whole galaxy. This formula 



is useful for distant (unresolved) galaxies that lie entirely 
within the spectrograph slit, or within the field of an 
integral-field spectrograph. For nearby galaxies, however, 
this formula is not useful since velocity-dispersion profiles 
typically go out only to the effective radius Re (eq. I14p . 
which contains only half the light. 



3.2 Circular speed from velocity dispersion for 
Sersic models 

3-2-1 Constant anisotropy 

Let us continue to assume that we have a spherical galaxy 
in a logarithmic potential and consider three variants of the 
stellar orbital anisotropy: isotropic orbits (i.e., the velocity- 
dispersion tensor is isotropic), circular orbits, and radial or- 
bits. It is straightforward to show from the Jeans equation 
(see Appendix [X| that the line-of-sight velocity dispersion 
in these three models is given by: 



a\R) = aLiR) = ' 
for isotropic orbits. 



i R 
■I{R) 



CT^{R) = Cr^irc(-R) = 2 L""^ 

for circular orbits and 



I{x) 



^LiR)] 



(11) 



(12) 



a (R) = crrad(-R) = -V, 



Rl(R) Jr 

for pure radial orbit^. In Figure [3] the radial dependence of 
the line-of-sight velocity dispersion is shown for the above 
three cases, assuming that the surface brightness of the 
galaxy is described by a Sersic law I{R) oc exp(— a7?^^") 
with index n — 1,4,12. The radius is plotted in units of the 
effective radius Re, the radius containing half the light: 



1 f° 
RI{R)dR = - / 

2 Jo 



RI(R)dR- 



(14) 



From Figure|3]it is clear that for all n and for all three models 
spanning the range of possible anisotropies the dispersions 
are rather similar at about 0.5-0.67?e. This result suggests 
that somewhere near this radius the line-of-sight velocity 
dispersion may provide a measure of the circular speed of 
the underlying isothermal potential that is relatively inde- 
pendent of the details of the stellar velocity distribution. 

There are many possible ways to define the "optimal" 
radius or "sweet spot" at which the relation between the 
line-of-sight velocity dispersion and the circular speed is 
likely to exhibit the smallest variation. In Figure |4] we plot 
the radii where the line-of-sight dispersions from each pair 
of the three models coincide, as a function of Sersic index n. 
From equations (|ll|l . (|12[) and (|13|) these are the solutions 
of the equations 

^ Here we ignore the possibility that some of the models con- 
sidered may have unphysical (negative) distribution functions at 
some radii. This is, for example, the case for a pure radial-orbit 
model when r — > so long as the emissivity di verges more slowly 
than j(r) oc jRichstone fc Tremaine|[l983) . 
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Table 3. Relative changes in circular speed with respect to the reference value x given in Table [2] when changes are made to the 
analysis procedure (see il2.4| l. 



Galaxy 


Aataund 


Err. 


Almxb 


Err. 


Ans 


Err. 


An x2 


Err. 




Err. 




/o 




% 




% 




% 




% 




ngcl399 


4.23 


0.80 


-2.01 


0.44 


2.36 


0.64 


1.40 


0.94 


1.22 


0.39 


ngcl407 


0.76 


5.67 


-1.61 


2.37 


8.25 


2.88 


5.66 


1.74 


1.05 


2.22 


ngc4472 


2.43 


2.17 


-0.98 


0.79 


-2.31 


1.16 


-0.42 


0.76 


-3.71 


0.70 


ngc4486 


3.88 


0.40 


-5.69 


0.37 


-4.58 


0.39 


1.44 


0.26 


-1.00 


0.31 


ngc4649 


-2.48 


1.14 


0.36 


0.41 


2.67 


0.62 


3.30 


0.43 


0.08 


0.39 


ngc5846 


1.53 


2.61 


-1.82 


1.02 


-2.57 


1.50 


-1.68 


0.81 


-5.14 


1.21 


Mean 


1.73 




-1.96 




0.64 




1.62 




-1.25 




RMS 


2.82 




2.69 




4.35 




2.89 




2.70 





Name 


changes in the analysis 


Section 


Aabuiid 


free metal abundance 




Almxb 


a power law is added 


(El 


Ans 


difference between North and South 




Arl 


ri X 2 


m 


Ar2 


r-2/2 


m 




Figure 3. Lino-of-sight velocity dispersion as a function of radius for spherical galaxy models with Scrsic index n = 1,4, 12 - left, 
middle, right panels respectively. The solid, dotted and dashed lines correspond to systems composed of isotropic, circular and radial 
stellar orbits. The thin vertical line shows the position of the "optimal radius" or "sweet spot" at which the deviations among the three 
models are minimized. 



nT{n, z) = z"e^^ circular = radial 

3z"nT{—n, z) = e^^ isotropic = circular (15) 

32^"'r(— n, z) = r(?i, z) isotropic = radial. 



where z = aRl''"', and F is the upper incomplete gamma 
function. Models of galaxy formation usually suggest that 
the velocity-dispersion tensor is isotropic near the center and 
radially biased (cr^ > ag) in the outer parts of the galaxy, 
so we should prefer the bottom half of the range of Rs indi- 
cated by these curves. A natural choice for the sweet spot Rs 
is the middle curve, where the line-of-sight velocity disper- 
sions for radial and circular models are equal (thin vertical 
lines in Figure [3]). This can be written for arbitrary surface- 



brightness profiles as 

poo 

/ I{x)dx = RsI{Rs). (16) 

J Rs 

The values of Rs, aiBo{Rs), o-circ{Rs) and a^adiRs) (the latter 
two quantities coincide by the definition of Rs) are shown 
in Figure [S] and in Table [J] as a function of Rs/Re- 

The range of dispersions spanned by these three mod- 
els is of course smaller than the range spanned by all possi- 
ble equilibrium models (see for example Figures 1 and 2 of 
iRichstone fc Tremaine|[l984l '). Our results nevertheless sug- 
gest that the line-of-sight velocity dispersion evaluated at 
Rs ~ 0.5i?e is a better proxy for the circular speed in a 
galaxy with a logarithmic potential than the dispersion at 
any other radius (e.g., the central velocity dispersion). 
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Figure 4. Radius Rsin) at which the line-of-sight dispersions 
(eqs. 111! 1121 and for systems composed from different types 
of orbits coincide, as a function of Sersic index n. Solid line; 
<^circ{R) = o"rad(^?), dashed line: criao(R) = a-rad(^2), dotted line; 

O-iBo(-R) = 0"circ(^?)- 



0.7 



0.65 




0.55 



Figure 5. The line-of-sight dispersions evaluated at Rs{n) [where 
o'circCJ^s) = frad(^s)] ^rc shown with solid (for systems with in- 
trinsically isotropic orbits) and dashed (radial and circular orbits) 
lines. 



Determination of the effective radius and the Sersic in- 
dex of real galaxies depends (sometimes strongly) on the 
range of radii used to fit the surface brightness profile and on 
the method of extrapolating to calculate the total flux. This 
is especially true for galaxies with extended envelopes, char- 
acteristic of central galaxies in clusters and groups. For such 
galaxies the reported values of Re and n for the same object 
may vary strongly. A strikin g example is M87: the effective 
radiu s of 145" in Table [1] ( D'Onofrio. Capaccioli. fc CaonI 
[l99i) is a factor of almost five smaller than the most re- 
cent value Re = 704" in iKormendv et al.l (|2009l ). Sim- 
ilarly t he reported values of the Sers i c inde x vary from 
6.51 in iD'Onofrio CapaccioU. fc Caoiil (|l994l ) to 11.84 in 
iKormendv et aL ( 20091 ) . The position of the sweet spot 
changes accordingly. Clearly such large variations have 
two main causes: (i) deviations of the observed surface- 
brightness profile from a pure Sersic profile, and (ii) the 
weak dependence of the surface-brightness slope on R/Re 
and n or, more precisely, strong covariance between Re and 
n when fitting to data. Indeed fo r the Sersic profile the slope 
is (e.g.. lGraham fc Driveij|2005h 



dlogljR) ^ / E 
dlogR ~ \Re 



l/n 



(17) 



Obviously when n is large the variations of slope at a given 
R are not very sensitive to the variations of Re- 

This raises the question of which set of n and Re to use 
and how strongly the final result (circular speed evaluated 
near the sweet spot) depends on these parameters. The def- 
inition of the sweet spot (eq. I16|) does not explicitly use the 
value of Re, but is instead most sensitive to the slope of the 
surface brightness profile; thus, for example, if I(R) oc R~^ 
for R > Ro then any radius R > Ro can equally well serve 
as the sweet point (see also i\3.3} . This implies that the eval- 
uated circular speed should not be sensitive to the precise 
values of Re and n so long as they provide a good fit to the 
surface-brightness profile and our other assumptions (e.g., 
circular speed independent of radius) are satisfied. 

To illustrate this point, let us assume that M87 has 
circular speed independent of radius, isotropic velocity- 
dispersion tensor, and surface-brightness profile that is ac- 
curately fit by a Sersic profile with the Kormendy et al. 
parameters, Re — 704" and n = 11.84, and that we esti- 
mate the circular speed using Table |3] and the D'Onofrio et 
al. parameters Re = 145" and n = 6.51. The resulting error 
in Ve will be only 3% even though the values of Re differ by 
a factor of almost five. 

In practice the sweet spot at large radius is less favor- 
able since measurements of the line-of-sight velocity disper- 
sion have only a limited radial extent. In particular, to esti- 
mate the circular speed using the effective radius from Ko- 
rmendy et al. one needs the stellar velocity dispersion at a 
radius of ~ 0.5i?e ~ 5-6', which is not known. For this rea- 
son we decided to use the effective Re = 145" of D'Onofrio 
et al. 

It is possible that equation [TB] does not have roots at 
radii where good measurements of the line-of-sight velocity 
dispersion are available. In this case one can look for the 
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Table 4. Dependence of Rs and line-of-sight velocity dispersions 
at Rb on the Sersic index n. Here Rs is such that a-circiR) = 
frad(K); Rs is measured in units of the effective radius Re, and 



cr is 


measured in 


units of the circular speed 


Ve. 


n 


Rs /Re 


O"iso 




(^iso ~ ^circ)/'^iso 


1 


U.09o8z4 


U.UOOOO I 


n ^Af\f\^'^ 


n 1 yi nc; un 


2 


n Q 1 


U.UIOUIU 


U.OO ( U\J 1 


n nnu/inriQ 
U.Uyo4UUo 


3 


U.ooUy47 


U.OUOD^U 


n Kaof,70 


U.U7U94bo 


4 


U.oluozo 






U.UO ( OOOD 


5 


u.oUoyio 




n c;R7c;Qe 
U. OD ( OOO 


r\ r\ A o o c A 

U.U4o5od4 


6 


n /I none; 1 


0.593877 


0.568907 


n f\A on A HQ 


7 


U. 491859 


0.591909 


0.569932 


0.0371287 


8 


0.486838 


0.590369 


0.570730 


0.0332661 


9 


0.482674 


0.589130 


0.571369 


0.0301473 


10 


0.479156 


0.588111 


0.571894 


0.0275738 


11 


0.476140 


0.587256 


0.572333 


0.0254124 


12 


0.473523 


0.586530 


0.572705 


0.0235703 


13 


0.471228 


0.585904 


0.573025 


0.0219810 


14 


0.469199 


0.585359 


0.573304 


0.0205953 


15 


0.467389 


0.584881 


0.573548 


0.0193760 


16 


0.465765 


0.584457 


0.573764 


0.0182946 


17 


0.464299 


0.584078 


0.573957 


0.0173289 


18 


0.462968 


0.583738 


0.574130 


0.0164609 


19 


0.461754 


0.583432 


0.574285 


0.0156766 


20 


0.460642 


0.583153 


0.574427 


0.0149643 


oo 


oo 


0.577350 


0.577350 






radius R with measured cr(i?) where the relative mismatch 
of the left and right sides of equation [T5] is smallest. 

In the limit of large Sersic index n, the surface bright- 
ness of a galaxy at ra dii not too far from the e ffective radius 
declines as (e.g. jGraham fc Driveijr2005f ) and therefore 
volume emissivity declines as r~^. While formally the total 
mass and the effective radius diverge as n goes to infinity, at 
the same time the line-of-sight velocity dispersion becomes 
inde pendent of the radius and the anisotropy of stellar or- 
bits (|Gerhard|[l993 ) as long as the shape of the velocity- 
dispersion tensor, ag/a^, does not depend on radius. In this 
limit cr^(_R) = const = u^/S at all radii (Table H]). 



3.2.2 Anisotropy changing with radius 

In the models considered above the anisotropy parameter 
was independent of radius. We now consider several sim- 
ple models in which anisotropy parameter /?(r) — 1 — 
ag{r)/ar{r) depends on radius according to 



r'^ + r^ ' 



(18) 



where /3i and /?2 are the anisotropy parameters at r = 
and r — >■ 00 respectively, is the anisotropy radius, and 
the exponent c controls the sharpness of the transition. The 
radial velocity dispersion is then 



cr^(r) 



jir)W{r) 
where 

W{x)^x^^'{x' + r:f^f'^-^'^^\ 



dx. 



(19) 



(20) 



A similar model with c = 1, /3i = 0.2, jS^ = 1 was used 
to des cribe the anisotropy profile in M87 bv lDohertv et all 
('2OO9I). Another special case of our parame trization is the 
Osipkov-Merritt model (|Osipkov|[l979l : [MerrittnlOSSi ). cor- 
responding to c = 2, Pi — 0, /32 = 1- 

In the left panel of Figure [6] we show the behavior of 
the anisotropy parameter /3 for c = 1, 2 and 4 and different 
values of the anisotropy radius ra. The Doherty et al. models 
(c = 1) are shown by black lines and the Osipkov-Merritt 
models (c = 2) by red lines, both for = 0.1-Re, IRe, 5Re. 
In the right panel of Figure[6]we show corresponding profiles 
of line-of-sight velocity dispersion for an n = 4 Sersic model. 
The legends of the lines are the same as in the left panel. For 
reference, the thick green lines show our standard isotropic, 
circular and radial models. 

We did not check explicitly that these solutions of 
the Jeans equation correspond to non-negative distribu- 
tion functions. Sersic models with P — 1, corresponding 
to purely radial orbits, certainly have negative distribution 
functions and thus are unphysical. For the Osipkov-Merritt 
models, the smallest anisotropy radius shown in Figure |6] is 
ra = 0.1i?e, close to but outside the range rg < 0.05-Re at 
which the distribution function is negative (jMerrittlllQSSh . 

Considering the models shown in the right panel of Fig- 
ure[6]it is clear that the "sweet spot" at which the sensitivity 
of the line-of-sight dispersion to the orbit structure is mini- 
mized is near 0.5-Re, just as we found in ii3.2.1l from examin- 
ing models with radial, circular, and isotropic orbits. How- 
ever, and not surprisingly, the models with radially varying 
anisotropy span a larger range in line-of-sight dispersion at 
the sweet spot, with the largest deviations occurring when 
ra/Re ~ 1. For instance at i? = 0.516-Re (taken from Table 
|4]for n — 4) the Osipkov-Merritt model with ra = Re pre- 
dicts a line-of-sight velocity dispersion 14% larger than the 
isotropic model and 20% larger than a model composed of 
circular or radial orbits. Making the transition in /3 smoother 
(c = 1) or sharper (c > 2) makes the wiggle in the line-of- 
sight velocity profile less (or more) pronounced as shown in 
Figure [6] 

So far we have discussed only sy stems with r adial 
anisotropy ("Type I" in the notation of iMerrittHlQSSl l. We 
have not discussed "Type 11" systems in which the orbits 
become predominantly circular at large radii, because such 
systems are not produced in our current model of galaxy 
formation. At least for the Osipkov-Merritt models (c = 2) 
these curves have an abrupt change in the l ine-of-sight ve - 
locity dispersion at R = ra (see Figure 5b in lMerritt|[l985l l. 



These curves fill the area between the isotropic and circular 
models and near R = Rs (n) they predict smaller dispersion 
than the isotropic model. 

Even though spatially varying anisotropy can increase 
the spread of the expected value of the line-of-sight velocity 
dispersion compared to our set of isotropic, radial and circu- 
lar models, the dispersion at the radius R ~ 0.5 -Re remains 
relatively insensitive to the orbital anisotropy and hence its 
value at this radius provides a useful measure of the circular 
speed of the underlying isothermal potential. 
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Figure 6. Left: Radial dependence of the anisotropy parameter /3. The Doherty et al. models (c = 1) are shown by black lines, and the 
Osipkov-Merritt models (c = 2) by red lines. The curves corresponding to different values of the anisotropy parameter ra (see eq. I18I I 
are shown with dotted, solid and dashed lines for ra = O.liJe, l^?e and 5-Re respectively. Right: Corresponding profiles of line-of-sight 
velocity dispersion for an n = 4 Sersic model. The legends of the lines are the same as in the left panel. The thick green lines show for 
reference our standard isotropic, circular and radial models. 



3.2.3 Aperture dispersions and other extensions 

The methods we have described for estimating the circu- 
lar speed Vc rely on the value of the line-of-sight velocity 
dispersion a-{R) at a single radius Rs. Obviously, one could 
imagine improved estimators based on some combination 
of the dispersions at two or more radii Rai,Rs2, .... A re- 
lated possibility is to use weighted averages of the form 
I{R)a^{R)f{R)dR/ I{R)f{R)dR for some suitably 
chosen function f{R). Here we discuss only one possibility, 
the use of aperture dispersions {f{R) = R), which we shall 
find to be less useful than dispersions measured at a single 
radius. 

Similar to the derivations in Appendix |^ one can 
find analytic formulae for the luminosity-weighted dis- 
persions inside an aperture of radius A, o"ap(^) = 
I{R)a^{R)RdR/ I{R)RdR for systems composed of 
isotropic, radial or circular orbits in a logarithmic potential. 
Thus 

, 1 , I(x)/x^dx 



circ(^) — ■^■'^c 



3 I{x)xdx 

, A'J^I{x)/x^dx 
6 I{x)xdx 



(21) 



2 , 1 2 , 2^/r-''(^)^^ 



,A^ I{x)/x^dx 



6 I{x)xdx 

In Figure [7] the dependence of the luminosity- weighted dis- 



persions on the aperture size A is shown for Sersic laws with 
index n — 1,4,12. The dispersions converge to the same 
value Wc/\/3 at very large aperture (as required by eg. [TO)l. 
but strongly diverge at small radii. For example for an n = 4 
Sersic model, the difference between aperture dispersions 
(crap,rad " cTap.circ) at A/R^ = 0.5, 1, 2 and 5 is 0.27, 0.17, 
0.10 and 0.034; for comparison the difference between the 
line-of-sight velocity dispersions aiso{R) — o'iad{R) evaluated 
at R = Rs{n = 4) is 0.034u<;. 

We conclude that the aperture dispersion, as a measure 
of the circular speed, is more sensitive to the anisotropy pro- 
file than the dispersion at a single "sweet spot" radius Rs 
unless the aperture is > 5Re. For spatially resolved galaxies 
this implies that the dispersion at the sweet spot is a more 
powerful tool for estimating the circular speed than the aper- 
ture dispersion. Nevertheless the aperture dispersion is often 
the only available quan tity which can be used to estimate 
the circular speed (e.g.. [Padmanabhan et al1l2004h . 



3.3 Circular speed from velocity dispersion using 
local measurements 

The methods outlined in the previous subsection require 
determining the galaxy's effective radius and Sersic index, 
which is often difficult, especially for the most massive and 
extended galaxies. Moreover, some galaxies are not well de- 
scribed by a Sersic law. In some cases one can use equations 
ITT}, (|12p or (|13p to predict the line-of-sight velocity dis- 
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Figure 7. Velocity dispersion within an aperture as a function of aperture radius for Scrsic index n = 1,4, 12. The solid, dotted, and 
dashed lines correspond to systems composed from isotropic, circular and radial stellar orbits. 



persion in terms of Vc without first fitting to a Sersic law; 
however, the use of these equations formally requires the 
knowledge of I{R) up to i? — > cxj. 

Let us now assume that we have an observed surface- 
brightness profile I{R), perhaps available only over a limited 
range of radii, and we wish to estimate the circular speed 
Vc from measurements of the line-of-sight dispersion cr(iZ) 
assuming that the potential is isothermal. One can differ- 
entiate equations (|11|) . (|12|) . or (|13p with respect to R to 
obtain a relation between the circular speed and the local 
properties of I{R) and a{R). Thus 



(^iBo{R) 



O-circ(-R) = 



1 -f a + 7 
1 a 



O-rad(-R) 



' 2 1-1- a -I- 7 
, 1 



(22) 



where 



dlogljR) 
dlogR ' 



7 ^ 



d log a 
' dlogR' 



d^ log[I{R)a^ 
d(logi?)2 



.(23) 



Although the terms 7 and 5 can be evaluated accu- 
rately if there is sufflciently good data on the dispersion 
curve ct{R), in a typical galaxy we expect these to be sub- 
dominant compared to a and . This is illustrated in Figure 
[8] for a galaxy with Sersic index n = 4. If one neglects the 
subdominant terms then the expressions simplify further to 



o-iso(-R) = Vc 



a + 1 



2 / D^ 2 1 



<^Ld{R) = «c 



2 a-f- 1 
1 



(24) 



For pure radial orbits we can also use the expression which 
neglects all derivatives of the line-of-sight velocity disper- 
sion, but keeps the second derivative of the surface bright- 
ness. Namely, 
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Figure 8. Various terms contributing to the relation between the 
circular speed Vc and the line-of-sight velocity dispersion, (see eq. 
1231 for definitions) . Each term is evaluated for an n = 4 Sersic 
model. The 7 terms are calculated for isotropic, circular, and 
radial orbits; the 5 term is calculated for radial orbits. The term 
a related to the slope of the surface-brightness profile (the upper 
curve) dominates. 



+ Si 



1' 



(25) 



where 5i = d'^ \og[I{R)]/d{\ogRf . 

In the case of a power-law surface-brightness law I{R) oc 
R~°' in which the anisotropy parameter /3 is independent of 
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Figure 9. Relation between the line-of-sight velocity dispersion 
and the circular speed for the case of a pure power law surface 
brightness distribution I{R) oc R~". 



radius the expression is (see also Figure [9]): 



1-/3- 



(26) 



a + 1-2/3 V " 1 + a, 

As pointed out by I Gerhard! (|l993f l. for a = 2 the anisotropy 
parameter cancels out. In this case the relation between cir- 
cular speed Vc and the line-of-sight velocity dispersion at 
a given radius R solely depends on the local slope of the 
surface-brightness profile (for (3 = const). 

The above considerations suggest that one can get a rea- 
sonable estimate of the circular speed using the observed sur- 
face brightness I{R) and line-of-sight dispersion cr(_R) pro- 
files over a limited range of radii, together with equations 
(|22p . Similar to equations H15|) we can introduce different 
types of sweet spots (hereafter R2) where pairs of orbital 
models with observed quantities cr{R), a, 7 and 5 corre- 
spond to the same circular speec|f|. In particular 



R2,ic '■ a = 2 isotropic = circular 

R2,ir '■ Of + 7 = 2 (if 5 = 0) isotropic = radial. 



(27) 



The expression for R2.cr (circular = radial) is more compli- 
cated (see eQ. I22|) but based on the example of a Sersic profile 
we expect that R2,ir < R2,cr < R2,ic so the range between 
the two radii in equation (|27p should include R2,cr- When 
both 6 and 7 are small then all three models intersect at 
the same point R2 where a = 2 and criso(i?2) = crcirc(i?2) = 

0-rad(-R2) = 

* Note that these sweet spots are not the same as defined in 



We illustrate the application of this procedure to 
two well-studied systems — M87 and NGC3379 (the first 
is in our X-ray sample but the second is not). Shown 
in Figures [10] and [11] are the observed surface brightness 
and line-of-sight dispersion profiles (panels A and B) to- 
gether with smooth curves approximating their large-scale 
trends. For M87, the photometry is from iKormendv et al] 
(j2009l ). The line-of-sight velocit y dispersion is a co mbi- 
nation of van der Marel data (Ivan der Marell [l994() for 



the central regions, SAURON (lEmsellem et al . 
the mid-radi al regions and VIRUS-P (iHill et al 



20041) for 



20081 ) 



presented in iMurhphv fc Gebhardtl (|2009|) f or th e outer- 
most stellar kinematics. Gebhardt fc Thomas! (|2009l ) present 
an analysis using the SAURON data, van der Marel 
data, and kinematics from globular cluster kinemati c s. For 
NGC3379 the photometry is from ICapaccioh et all (| 19901 ) 
and Gcjj hardt ct al. ( 2 00d). and th e vel ocity dispersion is 
fromlstatler fc Smecker^Hand (| 19991 ) and lKronawitter et al.l 
(|200d ). We also use planetary nebula velocity dispersions 
at large distances where stellar velo city dispersions are not 
avai lable. For M87 we use data from lArnaboldi et al] (|2004l ) 
and iDohertv et all (l2009h. and for NGC3379 we use disper- 
sions as calculated b y de Lorenzi et al.l (|2009h from data in 
iDouglas et al] (|2007l ). The use of these data in our analysis 
is based on the empirical result that the spatial distribution 
of pl anetary nebulae trac ks the spatial distribution of the 
stars ICoccato et al.l (|2009l ). but see the discussion below. 

We then use the interpolated curves I{R) and cr(i?) to 
calculate the derivatives a, 7 and 5 defined in equation (|23p . 
The values of a, 7, 5 and the sum a + 7 are shown in panel 
C with red, blue, green and black lines respectively. The 
term S, which depends on the second derivative of the sur- 
face brightness and the line-of-sight velocity dispersion, is of 
course most sensitive to the way the data are interpolated, 
but usually makes only a small contribution to the estimate 
of circular speed. Some of the wiggles in S seen in panel C 
of Figures [TOl and fTTI are certainly spurious. 

As we argued before there are two robust methods to 
estimate the circular speed (for a logarithmic potential): (i) 
fit the photometry to a Sersic or other model and determine 
the circular speed from the line-of-sight velocity at the sweet 
spot Rs f H3.4.ip : (ii) determine the circular speed from the 
range of radii between points where a ~ 2 and 7 -I- a « 2 
f i)3.4.2l) . At these radii the anisotropy of orbits does not af- 
fect the estimate of the circular speed (all three black curves 
intersect in the vicinity of this point). 

M87 is fit well by a Sersic model with n = 6.5, and the 
first method yields Rs = 72 arcsec and Vc — 1.7a{Ra) = 
530 km s^^ (Tables [1] and [4|. The second approach, based 
on radii ~ 125 arcsec, yields Vc ~ 536 km s~^ in good agree- 
ment. In principle one can extend this analysis to larger 
radii using planetary nebulae (PNe), the data points shown 
as red open symbols in Figures[TO]and[TT] However, there are 
two problems in determining circular velocities from these 
outer data points: (i) the PNe dispersions are reported for 
R ^ 800", leaving considerable freedom in the interpola- 
tion of a{R) between ~ 200" and 800"; (ii) at these large 
distances, the logarithmic surface brightness gradient a ap- 
plicable to the outer halo of M87 has some uncertainty, due 
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to both the uncertain details of the inferred truncation of 
the halo lDohertv et alj (|2009l ') and the possible contribution 
from the intracluster light (ICL) to the outer surface bright- 
ness profile. 

NGC3379 is fit by a de Vaucouleurs model (a Sersic 
model with n = 4) with Re — 50 arcsec (|Capaccioli et al] 
I1990I ). so the first method yields Rg ~ 26 arcsec and 
Vc — 1.7a{Rs) — 28f km s~^. The second approach, based 
on radii near 50 arcsec, yields Vc — 272 km s^^ in good 
agreement. For NGC3379 the stellar and PNe dispersions 
overla p in radius, and no ICL has b een found in the Leo 
group [Castro-Rodn'iguez et al.l (|2003f ). so we used the stellar 
surface brightness profile and the combined stellar plus PNe 
line-of-sight velocity dispersion profile (solid line in panel B 
of Figure [TT|) for the analysis. 

In the bottom panel (D) of Figures [TOl and [TT] we show 
the estimated values of the circular speed as a function of ra- 
dius, for various assumptions about the velocity anisotropy. 
Red points are derived from the values of the line-of-sight 
velocity dispersion, shown in panel B, which were converted 
to the circular speed using equation (|22|) for isotropic orbits. 



Vc{R) = rJopt{R) X ^/m^+7. 



(28) 



The black solid, dotted, and dashed curves are the result of 
applying equations (|22|l for isotropic, circular and radial or- 
bits respectively to the smoothed curve approximating the 
line-of-sight dispersion profile cr(i?) (black curve in panel B, 
see also Fig. lI2p . The agreement of the black solid curve with 
the red points in panel D is a measure of the errors intro- 
duced by smoothing. Our starting assumption was that the 
galaxy potential is logarithmic at all radii (i.e., the circular 
speed is constant), so the results in Figures [Til] and [TT] are 
self-consistent only if the circular speed shown in panel D is 
approximately independent of radius. In M87 the solid black 
curve is rising at radii larger than a few tens of arcseconds, 
while in NGC3379 the black curve is instead declining with 
radius. This implies that either the assumption of isotropy 
is incorrect — in which case the potential may still be log- 
arithmic, and the methods of iJSjwill still give an accurate 
assessment of the circular speed — or the assumption of con- 
stant circular speed is incorrect. As discussed in H3.3.Il if the 
orbits are indeed isotropic then the black curves in panels D 
of Figures (TDj and [TTj are approximately correct. In case the 
anisotropy of orbits is not known the estimates of the local 
circular speed are still robust, but only in the vicinity of the 
radius where — dlog7(i?)/dlog_R « 2. 

Thin blue lines in Figure [T^j show the circular 
spe eds from published detailed dyna mical analyses of M87 
by iRomanowskv fc Kochanekl l200ll . their models NFWI, 
NFW2 and NFW3. The green soh d line in Figure [T?l 
shows the best fitting model of ICebhardt fc Thomas! 
ll2009l). The difference in the circ ular speeds between 
IRomanowskv fc Kochane^ I2OOII and iGebhardt fc Thomas! 
l|2009t ) can be traced largely to the difference in the 
line-of-sight velocity dispersion data and to the absence 
of a central black hole in the IRomanowskv fc Kochane^ 
!200ll model. Our estimates of the circular speed, based 
on SAURON and VIRUS-P data of IGebhardt fc Thomas! 



(|2009l ): iMurhphv fc Gebhardtl (|2009l ) and the assump- 
tion of isotropic o r bits, agree reasonably well with 
iGebhardt fc Thomas! (|2009|) between 20 and 200 arcsec. 
Compared to Romanowskv fc Kochane^ (|200fl ) our results 
(the solid black line) predict significantly larger circular 
speed over the entire range of radii. This is not surpris- 
ing, since new measuremen ts (jCebhardt fc Thomad [2OO9I : 
iMurhphv fc Gebhardtl |2009| ) give larger line-of-sight veloc- 
ity dispersions. They also disagree inside f arcsec but here 
the Romanowsky fc Kochanek models exhibit substantial 
orbit anisotropy so we would not expect agreement with our 
isotropic models. In addition, the Romanowsky fc Kochanek 
models do not include a central black hole, which makes a 
significant contribution to the circular speed inside 2-3 arc- 
sec. Finally, the thick gray curve in Figure[T2]is the circular 
speed derived from heavily smoothed X-ray data. For this 
curve we used the potential profile obtained from XMM- 
Newton data (see Fig. [T|, which was smoothed using a a 
Gaussian filter similar to the filter described in eq. IBf I of 
Appendix [B] with Ai? = 0.5. We differentiate smoothed po- 
tential to obtain the circular sp eed — r dtp/ dr. As dis- 
cussed in lChurazov et aL (I2OO8I), the X-ray data agree we ll 



Romanowskv fc KochaneM (|200lf ) 



with the NFW2 model of 
This is also seen from the comparison of circular speeds in 
Figure [T2] - compare the thick gray line (X-rays) and the 
blue line, corresponding to NFW2 model (the curve with 
the largest circular speed at ~ fOOO"). The new optical data 
suggest substantially higher circular speed and a larger con- 
tribution of the non-thermal pressure is needed to bring the 
X-ray and optical data in agreement. We also note that all 
methods suggest the increase of the circular speed in M87 
from ~400-500km s"^ inside central 2' to ~600-700km s"^ 
outside fO'. 

In NGC3379 our prescription for evaluation of the cir- 
cular speed assuming isotropic orbits (black solid line) is 
remarkably consistent with detailed models (blue lines) be- 
tween radii of 3 arcsec and fOO arcsec, even through the 
circular speed in the detailed models is far from fiat. We 
note here that the models with the highest and lowest circu- 
lar speeds at 400" are formally ru led out by the likelihood 
analysis in lde Lorenzi et al] (|2009l ): the preferred models are 
the middle three. 

Alternatively one can neglect the contribution of 7 and 
5 and use equation (|24|l instead. The accuracy of the recov- 
ered value of Vc is illustrated in Figure (THj for several val- 
ues of the Sersic index n. We used this simplifie d approach 
to plo t the data for M87 globular clusters from lCote et al.l 
(|200lr ). We plot the line-of-sight velocity dispersions from 
Table f of Cote et al. in Figure [TO] with magenta squares 
(panel B). From Figure 12 of Cote et al., we estimated the 
power law slope of the surface density of globular clusters 
as a ~ f.35 and used equation (|24|l for isotropic, circular 
and radial models to evaluate Vc in panel D of Figure [10] As 
one can see the range of values of Vc spanned by these three 
models overlaps reasonably well with the results based on 
the stellar velocity dispersion, and the rise in Vc for r > 100 
arcsec appears to be mirrored in the magenta points. How- 
ever the rather shallow surface density of globular clusters 
I{R) oc R'^ "^^ leads to relatively large uncertainties in Vc- 
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Figure 10. Circular speed in M87 estimated from ti ie line-of-sigiit ve l ocity dispersion using equations I I22II . Shown in panel A is the 
observed surface-brightness profile in the V band from lKorinendj ;_et_al ] ll200gll (I(R) = WQ-i^^^, where ^ is in magnitudes per arcsec^). 
Shown in panel B is the line-of-sight dispersion profile from lDohertT^ al.l 1I2OO9I ). Black solid line s in pan els A and B show smooth fitting 
curves. The magenta points with large error bars in panel B are the globular cluster data from lCote ct al. ( 20 01]); these are shown for 
refe rence but are n o t used in any fits. The last two data points at 800-2000" are the PNe data from lArnaboldi et al.l 1I2OO4I ') (open circle) 
and lDohertv et al.l l|2009l ) (open triangle). The fitting curves were used to compute a = —dlog I{R)/dlog R, 7 = —d log cr^ (/?)/ci log i?, 
S = d? log[I(R)cr^]/d(log R)^ , and the sum a-l-7, which are shown in panel C by red, blue, green, and black lines. Shown in panel D are 
the estimated values of the circular speed. Red points in panel D correspond to the measured values of the line-of-sight velocity dispersion 
(shown in panel B), which were converted to the circular speed using equation I I22I1 for isotropic orbits and the fitted curve for a + -y 
in panel C. The black solid, dotted and dashed curves are the result of applying equations l|22|l for isotropic, circular and radial orbits 
respectively to the smoothed curves approximating the surface brightness I{R) and the line-of-sight dispersion a(R). The magenta points 
in panel D are the globular-cluster data converted to circular speed by approximating the surface density of globular clusters as a power 
law with a 1.35 and using equation I I24II for circular (stars), isotropic (solid squares) and radial (open squares) models. Thin vertical 
lines in panels B and D show the radii used to estimate the circular speed using the methods, described in il3.2.1l and i]3.3l respectively. 
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Figur e 11. The same as in F igure 1101 for NGC3379. Panel A sho ws the surface-b r ightne ss profile from wide-field B-band photome- 
trv by [Capaccioli et al ] lll99Cll) and HST V-band observation s from [Cebhardt et al.1 l|2000l V The velocity dispersion (panel B) is from 
IStatler fc Smecker-Ha!n3 lll999h and iKronawitter et all 1 I2OOOI) . supplemen ted by planetary nebu la velocity di spersions (points sh own as 
open squares with large error bars at /J > 30" in panel B, as calcul ated bv | de Lorenz i et al ] |2009l from data in lDouglas et a0l2007l ). Thin 
blue lines in panel D show the circular speeds from models A— E of lde Lorenzi et al. I boogll . 



3.3.1 Circular speed varying with radius 

Experiments in solving the Jeans equation for a galaxy with 
an n = 4 Sersic surface-brightness profile have shown that 
the local approximation works reasonably well at all radii 



even if the circular speed is not constant, provided that the 
orbital anisotropy is known and does not vary much with ra- 
dius. This is illustrated in the left panel of Figure [T4l Three 
types of model Vc{r) profiles are shown in this figure with 
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Figure 12. Circular speed in M87. Enlarged ve r sion o f panel D of Figure [TOl Blue lines show the circular speeds from the detailed 
dynamic al analysis of iRomanowskv fc Kochanekl l|200ll )— their models NFWl, NFW2, NFW3. The green line shows the best-fitting 
model of lGebhardt fc Thomaj ll2009h . The magenta points are the globular-cluster data converted to circular speed by approximating 
the surface density of globular clusters as a power law with a ft! 1.35 and using equation l|24|l for isotropic systems. The thick gray curve 
is the circular speed derived from heavily smoothed X-ray data. 



thick solid lines: (i) Vc{r) = const, (ii) Vc{r) steadily rising 
with radius and (iii) Vi.{r) having a maximum at r ~ i?^ 
and declining towards smaller and larger radii. For each of 
the Vc{r) profiles we calculated the projected velocity dis- 
persion ct{R), assuming constant anisotropy (/3 = 0,-9,1) 
and converted a{R) back to Vc{r) using eq. [22]for the ap- 
propriate type of anisotropy. The recovered Vc{r) curves are 



shown in Figure [14] (left panel) with thin solid, dotted and 
dashed lines for isotropic, circular and radial orbits respec- 
tively. The vertical gray line shows the radius R2 where 
—dlog I{R) /d\og R = 2, which is close to the expected loca- 
tion of the region where the sensitivity of recovered Vc{r) to 
the orbital anisotropy is smallest. Clearly the circular speed 
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Figure 14. Left: Recovery of the non-constant circular speed profile from tfie the hne-of-sight velocity dispersion profiles using cq. 
1221 when the anisotropy parameter is constant and known. Three types of model Vc{r) profiles are shown with thick solid lines: (i) 
Vc{r) = const, (ii) Vc(r) steadily rising with radius and (iii) Vc{r) having a maximum at r ~ Re and declining towards smaller and larger 
radii. For each of the Vc{r) profiles the line-of-sight velocity dispersion ct{R) was calculated for = 0, -9, and 1 and converted back to 
Vc{r) using eq. II22I I for the isotropic (thin solid), circular (dotted) and radial (dashed) orbits. The vertical gray line shows the radius R2 
dlogI{R) 

where = 2. The circular speed is recovered well even far away from R2, if we know the anisotropy of orbit distribution. Right: 

dlogR 

Ratio of the recovered circular speed profile to the true circular speed profile when the true Vc(r) is non-constant and the (unknown) 
anisotropy parameter varies with radius. The vertical gray line shows the radius R2. Without prior knowledge of the orbital anisotropy 
the circular velocity profile can not be recovered accurately, except in the vicinity of R2 where the spread in the Vc values (including 
models with variable Vc and /3) is small. 



is recovered well even far away from R2, if we can guess the 
anisotropy of the orbit distribution correctly. 

If the anisotropy varies with radius the error introduced 
in the circular speed is similar in amplitude to that for the 
case of a logarithmic profile (Figure[6l). Near the sweet point, 
defined according to equation H27|) . the value of Vc is recov- 
ered to within ~10-20% even if the circular speed is varying 
with radius. This is further illustrated in the right panel 
of Figure 1141 In this figure we use the same three mod- 
els of Vc{r) profiles as in the left panel and five different 
models for the orbital anisotropy: the constant anisotropies 
/? = 0, —9, 1, a model with orbits changing with radius from 
isotropic to radial (c = 1, /3i = 0, /?2 = 1, ''a = 2 in eq. 
lisp and a model with orbits changing from isotropic to cir- 
cular (c = 1, /3i = 0, P2 = —9, ra = 2 in eq. [T5)) . For each 
combination of Vc{r) and /3(r) we calculate the projected 
velocity dispersion cr{R) and convert it back to the circular 
speed eq. [22] assuming constant anisotropy /3 = 0, — cxj, 1. 
This results in a total 3 x 5 x 3 = 45 recovered circular 
speed profiles. The ratio of the recovered profiles to the true 
circular speed profiles is shown in Figure 1141 Not surpris- 
ingly without any prior knowledge of the orbital anisotropy 
the circular speed profile can not be recovered accurately. 



However, in the vicinity of R2 the spread in the Vc values 
(including models with variable Vc and /3) is smallest: the 
ratio Vc/vc,true varies between 0.78 and 1.12 for the whole 
set of models. 

At the radii where a is far from 2 (either smaller or 
larger) our method does not guarantee the recovery of Vc 
without a firm prior on the anisotropy parameter. At these 
radii additional sources of information such as higher order 
velocity moments (/13, hi) and more elaborate methods (like 
Schwarzschild's method) are needed. One example of a loca- 
tion where our method fails is the inner few arcsec in M87 
where a black hole of a few times 10^ Mq dominates the 
mass. Because the surface brightness is much fiatter than 

in this region the allowed range for the circular speed 
Vc (for anisotropy jS between —00 and 1) is very wide (panel 
D of Fig. I10|l , easily large enough to accommodate a black 
hole with a mass ranging from to 10^'' Mq or even largeiQ. 

^ Note that for the innermost {R. 55 1") points in the line-of-s ight 
velocity dispersion data (which are from lvan der Mare]||l993) the 
seeing and the size of the slit can strongly affect the dispersion 
measurement. No attempts to correct for these effects were made 
in the present study 
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Figure 13. Circular speed estimated from the line-of-sight ve- 
locity dispersion using approximate equation II 24 I I for isotropic 
(solid), circular (dotted) and radial (dashed) orbits, instead of the 
exact equation H22|l . The former equation uses a pure power law 
approximation to the surface brightness profile and completely ig- 
nores the terms 7 and <5 (sec i]3.3ll . For radial orbits more accurate 
values of the circular speed arc obtained (dash-dotted lines) from 
another approximate equation II25I I. Colors denote Sersic models 
with n = 4, 8, 12, and 20. 



3.4 Summary of suggested procedures for Vc 
determination 

3.J1-.I For a galaxy with the surface-brightness profile well 
described by a Sersic law 

Use Table U (column 2) to find the "sweet spot radius" 
Ra/Re for a given n; use the line-of-sight velocity dispersion 
at this radius and convert it to the circular speed using col- 
umn 3 of TableUl we call this Vc.s (see Table[T]). Use column 
4 from the same table to crudely estimate the uncertainty 
in Vc. 



3.4-2 For a generic surface-brightness profile 

Make a suitable interpolation of the surface-brightness pro- 
file (see, e.g., Appendix[B| and the line-of-sight velocity dis- 
persion profile. Calculate the slopes of these profiles (quanti- 
ties a, 7, 5, as defined in i)3.3|) . Convert interpolated line-of- 
sight velocity dispersion into circular speed using equations 
H22[l. Locate the sweet spot where the three curves are as 
close as possible (approximately corresponding to a range 
between radii where a ~ 2 and a -I- 7 « 2) and use the cir- 
cular speed value predicted by the curves in the vicinity of 
this radius. We call this Vc.i. When the curves do not inter- 
sect within the radial range well covered by the data (or in 



which, at least, the extrapolation is reasonable) one can use 
the radius where the range of circular speeds spanned by 
isotropic, circular and radial models is minimal. The differ- 
ence between Vc predicted by the three curves can be used 
to crudely characterize the uncertainty. 



4 RESULTS AND DISCUSSION 

Guided by the above considerations we compiled velocity 
dispersion measurements from the literature for the galaxy 
sample in Table [T] 

• For the central velocity dispersion (ctc) we rely on the 
Hyperleda database. The corresponding circular speed was 
estimated as Vc,c = V^ffc, an approximation often used in 
the literature but with little or no theoretical justification 
(column 9 in Table [T}. 

• For the "sweet spot" estimates Rg we 
use the values of the effective radius and the 

Caon. Capaccioli. fc D'Qnofriol 



(Il994h 



Sersic ind ex from 

(119931). iD'Onofrio. CapaccioU. fc CaonI 

iMahdavi. Trentham. fc TuUvl (|2005l ). and ISpolaor et all 
(|2008h . The values of Rs were determined from the effective 
radius and Sersic index using Tabled T he line-of-sight ve- 
locity d ispersions at R s were t aken from Kronawitter et ahl 



(|2000l ). ISpolaor et al] (l2008h. iDohertv et al 



references therein), and lOe Bruvne et al.l (|200ll 
corresponding values a{Rs) are given in column 8 of the 
Table [1] In practice, for the selected sample Rs is always 
close to Q.5Re and the ratio aiso/vc — 0.6. For one case 
when kinematic data do not extend to Rs (NGC4472) we 
used the outermost data points to evaluate a. Given the 
uncertainties in R^, n and the kinematics data, the accuracy 
of a is unlikely to be better than ±5-10 km s~^. In column 
10 we estimate the circular speed using equation (|ll|l for 
isotropic orbits (see also the entries for (Jiso in Table |4]) . 

• For a second estimate of the circular speed we use 
the procedure described in i)3.3l that is, we locate the re- 
gion in which the logarithmic slopes a and a -\- j (eqs. 
I23p are as close as possible to 2, find the line-of-sight 
dispersion a at this radius, and then estimate the circu- 
lar speed using equation (|22|l for an isotropic model, as 
Vcj — a {1 + a + j)^'''^ . The values of Vcj for each galaxy in 
our sample are given in column 11 of Table [T] The most un- 
certain is the value of Vc.i for NGC4 4 72, w here the kinematic 
data of iBender. SagUa. fc GerhardI \l994) end at i? ~ 40", 
whil e the Sersic approxima tion to the surface brightness pro- 
file (|Kormendv et al]|2009l ) suggests that the intersection of 



curves occurs at 



200" 



Following the discussion in ij2]we compare the circular 
speed derived from the X-ray analysis (Table [2| with the 
three estimates of the circular speed from the optical data 
as shown in Figure [151 Assuming that Vc,opt = 77 x Vc,x we 
derive the following best fitting relations: 



Vc,c = 1.12 X Vc,x RMS = 0.07 
Vc,s = 1.14 X Vc,x RMS = 0.11 
Vc I = 1.10 X Vc,x RMS = 0.09, 



(29) 
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Figure 15. Relation between ^cXi tlie estimate of the circular 
speed from X-ray data, and three forms of Vc,opt, the estimate 
from the optical data. Red squares - iic,c, blue triangles - Vc,s, 
black circles - v^^i- The dashed line corresponds to Vcfipt = ^c,x, 
and solid red, blue and black lines to the best fitting relations be- 
tween t>c,x and Vc,c, Vc,s and lie ; respectively. Small black squares 
shows the values of v^^x obtained when changes are made to the 
analysis procedure (see i|2.4l l. These changes include: letting the 
abundance of heavy elements be a free parameter; adding a power 
law to account for potential contribution of LMXBs; varying lower 
or upper limits of the radial range used for v,. x calculations and 
separate calculations of the circular speed for the Northern and 
Southern parts of each galaxy. These different values of v^- x are 
plotted at fixed Vc^Opt = ^c,c for each object. The black error 
bars show the crude estimate of the uncertainty in v,. x discussed 
in ^2A\ 

where root-mean-square deviations were evaluated as 

RMS = [J2f^,{v^.opt,^/v^,x.i-v}ViN where 
summation is over the objects in the sample, N — 6. We 
note here that the smallest scatter in the above relation is 
found for the circular speed estimated from the central ve- 
locity dispersion Vc,c = \/2o"c, rather than for more sophisti- 
cated estimates, which were constructed to be least sensitive 
to unknown orbit anisotropy. However, given the small size 
of our sample, no firm conclusion on the quality of different 
proxies to the circular speed is possible. 

When five of the six objects are analyzed the value of 
?7 varies from 1.07 to 1.13 for Vc.i (compared to 1.10 for 
the full sample of six objects). The conclusion that 77 > 1 
continues to holds. One can also estimate the significance of 
the statement that 77 > 1 from the fact that almost all of the 
points are in Figure [T^ are in the upper left part of the plot. 
Assuming that 77 can equally probably be larger or smaller 
than unity the probability that 6 out of 6 points (for Vc,c) 



or 5 out of 6 (for Vc,3 and Vc.i) are found with ry > 1 is 0.016 
and 0.11 respectively. 

As discussed in H2. 41 changing the assumptions made in 
the analysis of X-ray data affects the value of Vc,x- This 
includes the effect of variable metalicity, contribution of 
LMXBs, dependence on the radial range used for evaluation 
of Vc,x and independent calculations of the circular speed 
for the Northern and Southern parts of each galaxy. The 
values of Vc,x obtained in each of the tests described in H2.4I 
are shown with small black squares in Figure [TS] For clarity 
these values of Vc,x are plotted at the same Vc,Opt = Vc,c- 
The black error bars show the crude estimate of the uncer- 
tainty in Vc,x of ~7% obtained in H2.4I under simplifying 
assumption that the contributions of different effects co-add 
quadratically. 

As discussed in COS the derived value of Vc,x is sensi- 
tive to non-thermal forms of pressure supporl[f|. This non- 
thermal support can be parametrized by its impact on the 
relation between Vc,x and the true circular speed Vc. 

vl = vlx (1 + ^) ' (30) 

where Pn and Pg are the non-thermal and gas thermal pres- 
sures respectively. Thus the observed relation (eq. I29|l pro- 
vides a measure of Pn/Pg: 

77^ - 1 0.21 - 0.29, (31) 

Pg 

depending on the adopted estimator of the optical velocity 
dispersion in equation (|29p . 

There are a number of obvious caveats (both on the X- 
ray and the optical sides) associated with the above analysis. 
In particular there is considerable uncertainty in the heavy 
metal abundance determination, mentioned in ij2l which 
may affect the estimates of the potential at the level of a few 
per cent. Some of our galaxies sit in massive group/cluster 
halos; in such cases the circular speed increases at large radii 
(e.g. M87 - see Figure I12|l , so the approximation of the po- 
tential as a logarithmic law characterized by a single circular 
speed may not be sufficient for accurate comparison of the 
X-ray and optical data over several decades of radius. Minor 
deviations from hydrostatic equilibrium, visible as wiggles in 
the potential profiles, could also contribute to the scatter. 
Nevertheless we believe that our main result — that the depth 
of the potential well derived from X-ray data (from hot gas) 
is systematically shallower (by few tens of percents) than 
the corresponding optical value (from stars) — is robust. As 
mentioned in COS these values are consistent with the cur- 
rent paradigm of AGN controlled gas cooling/heating in the 
centers of clusters and indivi dual massive elliptical galax- 
ies (e.g.. IChurazov et al.ll200^ '). The key assumptions of this 
paradigm are that: 

• Gas radiative cooling losses are balanced by the me- 
chanical energy provided by a central black hole (AGN); 

• Dissipation of the mechanical energy is occurring on 

^ In COS instead of notation = Jf ^ used = 
2(T^ logr, where v^^x = V^crx- 
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time scales comparable (within a factor of few) to the sound 
crossing time of the cooling region; 

• The ratio of the dissipation and cooling time scales (of 
the order of 0.1-0.2) sets the ratio of the non-thermal and 
thermal energy densities. 

We finally note that the sample considered here is not 
statistically complete and generalization of these results to 
all gas-rich ellipticals must be done with caution. 



5 CONCLUSIONS 

Using non-parametric method we reconstructed gravitating 
potentials for a sample of six X-ray bright elliptical galaxies 
observed with Chandra and XMM-Newton. Their gravita- 
tional potentials can be reasonably well approximated by the 
isothermal (logarithmic) law Vc.x log r + b over the range of 
radii from ~0.5 to ~25 kpc, corresponding to ~ 0.05-3 Re- 
This result is in line with recent lensing data which also 
suggest isothermality of the gravitational potentials of the 
early-type gala:xies and earlier results based on stellar kine- 
matics and X-ray data. Many galaxies in our sample are 
located at centers of massive groups/clusters and the X-ray 
data going beyond the optical extent of galaxies show the 
steepening of the potential at large radii. 

We suggest two new methods to derive an estimate of 
a galaxy's circular speed if its potential is described by the 
isothermal law. (i) For a spherical galaxy with a given Sersic 
index n, the line-of-sight stellar velocity dispersion aopt{Rs) 
evaluated at Rs = -Rs('^) ~ 0.5-Re turns out to be rela- 
tively insensitive to the anisotropy of stellar orbits, (ii) We 
suggest the extension of this method for a generic surface- 
brightness profile which allows one to estimate the circular 
speed through the "local" values of the line-of-sight veloc- 
ity dispersion and the logarithmic derivatives of the surface 
brightness and the velocity-dispersion profiles. 

Application of these methods to a sample of six mas- 
sive elliptical galaxies and the comparison to the results of 
X-ray analysis suggests that on average 20% of the gas pres- 
sure in these galaxies is provided by non-thermal compo- 
nents (e.g., micro-turbulence or cosmic rays). This result 
is in broad agreement with the current paradigm of AGN 
controlled gas cooling/heating in the centers of clusters and 
individual massive elliptical galaxies. 
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APPENDIX A: LINE-OF-SIGHT VELOCITY 
DISPERSION IN A LOGARITHMIC 
POTENTIAL 

The goal of this Appendix is to derive formulae for the 
line-of-sight velocity dispersion profile <j{R) for a spherical 
galaxy with surface- brightness profile I{R), assuming that 
the gravitational potential is logarithmic. 



log r + const. 



(Al) 



The line-of-sight dispersion will depend on the shape of the 
velocity-dispersion tensor, defined by its radial and tangen- 
tial components (Jr(r) and (^'^{r) = o-g(r). We examine three 
simple cases that should span the range of possible behav- 
iors: (i) isotropic orbits {a^ = = (ii) radial orbits 
(cr^ = CTg = 0); (iii) circular orbits (ct^ — 0). 

In a spherical system, the volume emissivity j{r) and 
the surface brightness I{R) are related by 



I{R) = 2 



rdr 

1 f°° dR dl 



T^Jr VR^-r^ dR' 
If the system is isotropic, the Jeans equation reads 



(A2) 



d 



.d4> 



(A3) 



Specializing to the logarithmic potential and integrating 



j{r)at.{r) = vt I —j{u) 



(A4) 



The line-of-sight dispersion at projected radius 7?, (Jiso{R), 
is given by 



I{R)aL{R) = 2 



r dr 



2v, 
2v. 



2 1 r dr 



j{r)(jl{r) 



du 



2 ; du\/ V? — 

J R 



u 



(A5) 



Replacing j{u) from equation (|A2|l and exchanging the order 
of integration. 



^ ' ^ ' TV J J, dxjj^ u^X^ - 



(A6) 



The inner integral is ^nil — R/ x) so after integrating by 
parts 



I{R)oL{R) = Rvl 



I{x) 



dx, 



(A7) 



which is equation (|lip . 

If the orbits are circular, a set of stars with random 
orientation at radius r and projected radius R contributes 
a line-of-sight dispersion ^VcR^/r^. Thus the line-of-sight 
dispersion at R is given by 



J R 



dr 



(A8) 
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Replacing j(r) from equation HA2|) and exchanging the order 
of integration, 



I{R)a!UR) 



dr 



(A9) 



The inner integral is ^■k/{Rx) so after integrating by parts 



2/(i?)y« x^ 

1 2 1 2 / 



dx 



(AlO) 

which is equation (|12|) . 

Finally, if the orbits are radial, the Jeans equation reads 

_rja^ = -jr-. (All) 

Specializing to the logarithmic potential and integrating, 
»,2 roc 



j{r)ar (r) 



duuj{u). 



(A12) 



If the orbits are radial, a set of stars at radius r and projected 
radius 7? contributes a line-of-sight dispersion al{r){r^ — 
R'^)/r^. The line-of-sight dispersion at projected radius R, 
CTradiR), is then given by 



I(R)al^{R) = 2 / ±XlL^j^r)aXr) 



2vi 



2 f"" drVr2 - _R2 



duuj{u). (A13) 



Replacing j(it) from equation (|A2|I and exchanging the or- 
der of integration. 



2vl f°° drVr^ - R^ 



f dxsj x"^ — r^ — 

Jr dx 



2vl r°° drVr^ - R^ 



X dx 



2v^ 

— - I dxxl(x) 



dr \Jr''- — R? 



I{x) 



(A14) 



The inner integral is jit{x — R )/{Rx ), so 

,,2 rca 

I{R)oLd{R) = i^j^ I{x){l~R^/x^)dx, (A15) 
which is equation (jlSp . 



set of surface brightness measurements I{Ri),i = 1, . . . ,n. 
Given that we are interested primarily in logarithmic deriva- 
tives it makes sense to use fi = log Ri and Si = logJ(iii) 
instead of R and /. We can then fit the linear relation be- 
tween fi and Si as Si — afi + b, where a and b are the free 
parameters of the fit. We want to find "local" values of a 
and b in the vicinity of a given radius Tio, which are based 
on a set of measurements not far from _Ro- We do that by 
choosing a weight function 



W{Ro,R) = exp 



(logi?o-logi?)' 



(Bl) 



where the parameter A_r controls the width of the weight 
function. The best fitting parameters a and b (minimizing 
RMS deviation) for a given value of /o = logi?o are given 
by 

E f^w^Si E - E WiS^ E f^Wi 



where Wi = W{Ro, Ri), the sums are over i 

E,mSi-aE,/.m 



biRo) = 



(B2) 
1, . . . , n, and 
(B3) 



The smoothed function i{R) = g"'^' '"^^ ^+''(^) is shown in 
panel A in Figures [10] and [11] with a black solid line. We 
used A_R — 0.3 for both objects. 



APPENDIX B: SMOOTHING THE 
SURFACE-BRIGHTNESS PROFILE 

The procedure for evaluating the circular speed that is de- 
scribed in ij3]3] requires calculations of the derivatives of the 
surface brightness I{R) and the line-of-sight velocity disper- 
sion a{R). The derivative of the surface brightness is ex- 
pected to be the most important term (see Fig. [HJ. The 
observed surface-brightness profiles typically have enough 
data points with small errorbars to make the selection of 
a simple analytic model difficult. To counter this problem 
we used a simple procedure to smooth the data. Consider a 
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